Hole-Pairs in a Spin Liquid: Influence of Electrostatic Hole-Hole Repulsion 
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The stability of hole bound states in the t — J model including short-range Coulomb interactions is 
analyzed using computational techniques on ladders with up to 2 x 30 sites. For a nearest-neighbors 
(NN) hole-hole repulsion, the two-holes bound state is surprisingly robust and breaks only when the 
repulsion is several times the exchange J. At ~ 10% hole doping the pairs break only for a NN- 
repulsion as large as V ~ 4J. Pair-pair correlations remain robust in the regime of hole binding. The 
results support electronic hole-pairing mechanisms on ladders based on holes moving in spin-liquid 
backgrounds. Implications in two dimensions are also presented. The need for better estimations of 
the range and strength of the Coulomb interaction in copper-oxides is remarked. 

PACS numbers: 74.20.-z, 74.20.Mn, 75.25.Dw 
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The discovery of superconductivity in Sr\^ x Ca x 
CU24O41 0| under high pressure has triggered a consid- 
erable effort to understand the physics of copper-oxide 
ladder materials. Early theoretical studies predicted the 
existence of superconductivity based on a purely elec- 
tronic mechanism ||. The rationale for the hole-hole 
attraction is that each individual hole distorts the spin- 
liquid ground state of the undoped ladder, increasing lo- 
cally the energy density. The damage to the spin ar- 
rangement is, thus, minimized if holes share a common 
distortion. Several calculations, including computational 
studies for the t — J model [|j-f|, have confirmed that 
an effective attraction between holes exists in environ- 
ments with short-range AF correlations. The fact that 
the hole-pairs transform as d x 2_ y 2 under rotations, i.e. 
the same channel in which the real doped 2D cuprates 
superconduct, provides additional support to the elec- 
tronic mechanisms described above. 

However, studies of hole pairing in spin liquids are usu- 
ally performed without the introduction of an intersite 
hole-hole Coulomb repulsion. The accuracy of such an 
approximation has been much debated, but there are few 
robust calculations supporting the various points of view. 
The neglect of electrostatic interactions is particularly 
important for values of J/t such as 0.3-0.4, presumed to 
be realistic, since in this regime the pair size is small, 
roughly 2a (where a is the lattice spacing) Never- 
theless, note that in the real materials at least the 5 d- 
orbitals of Cu, as well as the 3 p-orbitals of O, should be 
taken into account. In this extended model, polarization 
effects caused by electron-hole excitations can effectively 
reduce the strength of the electrostatic repulsion. Actu- 
ally, 2 ~ has a large polarizability that influences on the 
effective Coulomb interactions of a variety of ionic insu- 
lators Q. This effect has been analyzed for the on-site 
repulsions Ud and U p [||. Although the bare couplings 
are substantially reduced by polarization, their strength 
remains the dominant scale and models with no doubly 



occupied sites capture properly this effect. However, it is 
unclear if the additional neglect of the hole- hole repulsion 
at distance a remains a good approximation. 

Unfortunately, the explicit calculation of the strength 
of screened Coulombic interactions is difficult. Using re- 
sults for 2D copper-oxides as guidance, a constrained- 
density-functional approach where the LDA bands are 
identified with a mean-field solution of the Hubbard 
model reports a repulsion U p d ~ 1.2eV between a pair 
of holes at the Cu and O ions ||. Since the bare value is 
7.8eU, a reduction of a factor 6.5 is effectively achieved 
by polarization effects. If the trend continues, then the 
repulsion V at distance a will be ~ 0.6eU ~ 4J (with 
J ~ 0.15eU (To) ) . However, other estimations compar- 
ing Auger spectroscopy results with (CU4O5) 6- cluster 
calculations using a multiband Hubbard model report 
Upd < leV [pf. In addition, for two holes in neigh- 



boring oxygens U p 



was observed 



suggesting 



that the hole-hole correlations decay rapidly with dis- 
tance. For this reason the above estimation 0.6eU should 
be considered as an upper bound for V. Other stud- 
ies based on two neighboring Zhang-Rice singlets report 
V ~ 0.2eV flljj . Finally, approaches where the dielec- 
tric constant e ~ 30 jl^] is used even at short distances 
provide V ~ 0.13eU. Then, current estimations locate 
Ve[ J, 4 J] as the realistic range for the hole-hole repulsion 
in a one band model. Such strong repulsion is dangerous 
for real-space AF theories. Coulomb interactions are also 
a central ingredient of the "striped" scenarios Q . 

The purpose of this paper is to discuss a computation- 
ally intensive calculation of the effect of intersite hole- 
hole repulsion on the hole bound states of the t — J model 
on ladders. Since the size of the hole pairs on planes 
and ladders is comparable for the same J/t, our results 
have consequences also for studies in 2D. The calculation 
is performed with the DMRG technique H], as well as 
using a diagonalization technique in a reduced Hilbcrt 
space [Jl5| (the Optimized Reduced-Basis Approximation 
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or ORBA) |ll|. The Hamiltonian employed here is the 
t — J model at a realistic coupling J/t — 0.4, supple- 
mented by a hole-hole repulsion Vhh = — j|) n i n j; 
where rii is the hole number operator at site i. The range 
R of the repulsion has been restricted to 1 and \/2 lattice 
spacings, since the speed of convergence of the numeri- 
cal techniques (both variational) decreases as R grows. 
Earlier results have been obtained in the case of static 
holes Q or using small clusters jl8| . 

Let us discuss first the case of two holes. Fig.l contains 
the hole-hole correlation C(j) = (noUj) obtained on 2 x N 
clusters with two holes and open boundary conditions 
(OBC) using DMRG (m = 300 states), for the case where 
Vhh acts only at a distance of one lattice spacing (i.e. 
R = 1) [lSf] . is a site at the center of the cluster, 
and the figure shows the hole-hole correlation along the 
leg opposite to where is located (results for the other 
leg are similar). C(j) is related with the probability of 
finding one hole at j when there is one at 0. 
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FIG. 1. Hole-hole correlation C(j) (see text) for 2 x 20 
(open circles) and 2 x 30 (full circles) clusters with two holes, 
using DMRG with OBC and J/t = 0.4. Results for several 
couplings V/J are shown. The Coulomb repulsion has range 
R = 1. Here j is the distance between a site at the center 
defined as origin, and sites belonging to the opposite leg. 

Figs.la-b correspond to V = 0, 2J and 4J. Here the 
results change only slightly as the lattice grows, and the 
two holes remain close to each other indicating the exis- 
tence of a bound state. Apparently a NN-repulsion al- 
ready larger than J is not enough to destroy the bound 
state, although it weakens it. On the other hand, Figs.lc- 
d show similar results but now for V = 8 J and 10 J where 
a substantial change in the hole distribution is observed 
as the cluster grows. The spreading of the hole over the 
entire lattice suggests either the breaking of the pair or 
a weak bound state. In the large-V regime one of the 
holes acts as a sharp "wall" to the other, which spreads 
its wave function in an effective square-well potential. 
Fig. 2a shows an average hole-hole distance defined as 



( d ) = Ej#o d j( n oWj)/Ej^o( n o"j)> where dj is the 0-j 
distance. The convergence of the results for V < 4 J as 
the size grows is clear. However, for V > 6 J, there is no 
obvious convergence. The curvature change of the 2 x 30 
results at V ~ 5 J also favors the interpretation that the 
bound state exists up to that coupling, since in the bulk 
(d) vs V/J will grow with positive curvature, diverging 
at the pair-breaking critical coupling. 
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FIG. 2. (a) Mean distance between holes (d), defined 
in the text, vs V/J obtained with the DMRG method at 
J/t — 0.4. Results for three lattice sizes are shown. The 
range of the Coulomb repulsion is R = 1; (b) Same as (a) 
but now using the ORBA technique with ~ 10 6 states, and a 
2 x 20 cluster with periodic boundary conditions. 

To further support these conclusions, the ORBA tech- 
nique has been implemented on a 2 x 20 cluster with 
periodic boundary conditions (PBC). The basis selected 
to optimize the convergence as the Hilbert space grows 
is the rung-basis, which contains 9 states for the t — J 
model [jl6|] . A convergence in the energy with up to 3 sig- 
nificant figures is achieved in this cluster by using ~ 10 6 
states in the subspace of zero momentum. While this is 
not as precise as the DMRG method for this model, the 
accuracy is enough for our purposes and ORBA has the 
advantage that it is performed with PBC. Fig. 2b shows 
(d) vs V/J. Here is any site of the lattice since there 
is translational invariance. For V / J < 4 it was observed 
that the hole-hole correlations change by a very small 
percentage when the space grows from 10 5 to 10 6 states, 
using as a starting configuration two holes in the same 
rung. In this regime the hole-pair bound state is tight. 
However, for V / J > 6 the changes in the correlations as 
the number of states grows are large. Actually the hole- 
hole distances reported in Fig. 2b for large V/J are not 
stabilized (thus they carry large error bars). Neverthe- 
less, as V/J increases the qualitative behavior of (d) is 
in agreement with the DMRG data. 

Let us now analyze what occurs when the range of 
the repulsion grows. Fig. 3 contains DMRG results for 
a repulsion of strength V and V/\/2 when holes are at 
distances 1 and \[2 lattice spacings, respectively, and zero 
otherwise. Although qualitatively similar to Fig.l, now 
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the region where the bound state is stable is restricted to 
V/J < 2, while results for V/J > 4 suggest the absence 
of a bound state due to the spreading of the second hole 
when one is at the center of the ladder. Fig. 4a shows (d) 
for the same interaction used in Fig. 3. In agreement with 
the previous discussion, the critical region is V / J ~ 3. 
For smaller (larger) couplings (d) converges (diverges) 
as the cluster size grows. ORBA results also suggest a 
critical coupling V/ J ~ 3. 
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FIG. 3. Same as Fig.l but with a longer ranged repulsion 
of strength V (V/\/2) at distance of one (\/2) lattice spacings. 

To complement the analysis using hole-hole correla- 
tions, Fig. 4b contains the binding energy defined as 
A B = E(2) + E(Q) - 2E(l), where E(n) is the lowest 
energy in the subspace of n holes. A negative im- 
plies that a bound state exists. At J/t = 0.4 and V = 
there is pairing of holes on 2-leg ladders |^,[| . The figure 
shows that for a NN-repulsion, the binding effect contin- 
ues even up to large values of V/ J. Actually the region 
6 < V / J < 10 corresponds to a weak bound state. The 
same figure shows results for range- V2 repulsion. In this 
case the critical coupling V/J\ c at which the bound state 
is lost is between 3 and 4. Both results are in agreement 
with the estimations based on Figs.3-4a. 

The previous analysis indicate that the stability of the 
two- hole bound states in t — J- like models depends on the 
value and range of the electrostatic Coulomb interaction 
between holes. When the Coulombic term is restricted to 
a realistic range R = 1, a repulsion as large as 0.6eV ~ 4 J 
weakens but does not destroy the pair, implying that the 
effective range of attraction caused by spin polarization 
is larger than one lattice spacing. Retardation effects 
(fully considered in the present calculation) due to the 
different energy scales of spin and charge excitations ( J 
vs t) likely contributes to the strong stability of the bound 
states Q. Note also that using, e.g., V/J = 4, the pair 
size is ~ 4— 5a (Fig. 2a). This result is close to estimations 
of the Cooper pair size in the 2D cuprates which report a 
coherence length £ c /a ~ 4 and ~ 3, for optimally doped 



La - 214 and YBCO, respectively ff(],|l| (using a = 
3.8A). To the extend that £ c /a are similar on planes 
and ladders, apparently a realistic NN hole-hole repulsion 
can actually improve quantitatively the predictions of the 
t — J model, without destroying the pairs. 
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FIG. 4. (a) Same as Fig. 2a but with a longer ranged 
repulsion of strength V (V/V2) at distance of one (v2) lattice 
spacings; (b) DMRG binding energy of two holes obtained for 
J/t — 0.4. Open circles and squares correspond to 2 x 20 and 
2 x 30 clusters, respectively, both with a Coulomb repulsion 
of range R = 1. Full circles and squares also correspond 
to 2 x 20 and 2 x 30 clusters, respectively, but now with a 
Coulomb repulsion of range R= y/2. 

The results discussed thus far have been restricted to 
the two holes sector. Since the pairs are very tight it 
is reasonable to assume that even with a finite (small) 
hole density the bound states will remain stable upon 
the addition of a NN-repulsion V. This assumption can 
be checked explicitly using a 2 x 20 cluster (shown to be 
already close to the bulk limit, see Fig. 4b). In Fig. 5a the 
binding energy defined as Ab = E(n)+E(n—2) — 2E(n— 
1) vs the hole concentration x = n/40 is shown for three 
values of V. As expected, the results are qualitatively 
similar to those of the two holes problem. While for 
V = the pairs remain stable even up to 20% doping 
and beyond, using a repulsion as large as V — 2J the 
stability persists roughly up to x ~ 0.17, and for V = 
4J up to x ~ 0.10. Finding a critical coupling V/J\ c 
decreasing with doping is reasonable since the stability 
of the pairs at finite x depends not only on the hole-hole 
attraction, but also on the mean distance between holes 
that introduces limitations on the pair size. As example 
consider V/J = 2. From Fig. 2a the pair size is (d) ~ 3. 
The density at which the hole mean-distance is also ~ 3 
corresponds to x ~ 0.17, in good agreement with the 
critical density at this V (Fig. 5a). Then, from the two- 
hole problem information about finite x can be obtained. 

Superconducting pair-pair correlations C(r) have 
also been measured here. C{r) is proportional to 
(A(m)AT(m + r)), where m, r are rung indices, and 
A(m) destroys a pair of electrons in a spin singlet at 
rung-m. An average over the whole cluster is used to 
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calculate C(r). At V — previous results were repro- 
duced as a test | f22[ . The new results are shown in Fig. 5b 
using hole-density x — 0.1 as example. The data is pre- 
sented at three characteristic distances (r = 3, 6, 10) |23f| 
vs V/J. The normalization using C(0) is important since 
the signal for superconducting correlations has an over- 
all penalization due to the reduction of the probability of 
finding two holes in the same rung as V/ J grows. This 
reduction only indicates that pairs become larger as V/ J 
increases. Fig. 5b shows that the strong V = supercon- 
ducting correlations remain robust as V/ J grows. Even 
for r as large as 10 the correlations are not negligible in 
the scale of the plot in the region < V/J < 4, which 
coincides with the region of hole binding (Fig. 5a). Sim- 
ilar conclusions were reached at x = 0.2 (not shown). 
Then, our results are compatible with a picture where 
hole-binding and strong superconducting correlations co- 
exist. Note, however, that as x — > the pair correlations 
must decrease due to the reduction in the density of pairs 
even if the binding energy is robust. 
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FIG. 5. (a) Binding energy of hole pairs obtained working 
at a finite density of holes on a 2 x 20 cluster. Results at 
J/t — 0.4 and several V's are shown; (b) Pairing correlations 
C(r)/C(0) (see definition in text) vs V/J at three character- 
istic distances r along the legs. J/t = 0.4 and the hole density 
is x = 0.1. The cluster size is 2 x 30. 

Summarizing, it has been shown that (1) the two- holes 
ladder bound state is more stable than expected upon 
the introduction of a NN hole repulsion; and that (2) 
this conclusion persists in the presence of a finite hole 
density. Hole pairs and robust superconducting correla- 
tions were found in the same region of parameter space. 
These results provide support to ladder theories that pre- 
dict hole-pairing based on electronic mechanisms that de- 
scribe holes as immersed in spin- liquid backgrounds 0. 
It also reinforces striped scenarios for cuprates where 
pairing is produced by carriers moving from the fluctuat- 
ing stripes to the ladders between them Regarding 
2D systems our results show that the effect of a Coulomb 
interaction is a subtle quantitative problem. Actually the 
binding energy of two holes on a 4- leg ladder with V = 
and J/t = 0.4 has also been estimated here. The result 
is |A^ L | ~ 0.14i (using a 4 x 12 cluster, m ~ 400 states, 



and truncation error ~ 10 5 ) which is similar to previous 
results in 2D (Fig.24 of Ref. {§). Then, the influence of 
the NN Coulomb interaction on planes will be stronger 
than on ladders. In this respect it is imperative to obtain 
either experimental or theoretical information about the 
range and strength of the Coulomb interaction in effec- 
tive one-band models for cuprates, specially in 2D. 
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